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ABSTRACT 

Infectious bursal disease virus (IBDV) is classified according to the antigenicity and virulence into classical 
virulent (cv), very virulent (w), and antigenic variant strains. The molecular basis for the IBDV antigenic 
variation is well established and is associated to the capsid protein, VP2 (gene VP2 of segment A), whereas 
both VP2 and the RNA-dependent RNA polymerase, VP1 (gene VP1 of segment B), have been correlated 
with the virulence. In this study, seventeen Brazilian IBDV samples previously characterized by the VP2 
gene as cv (three) and vv (fourteen) strains were genetically and molecularly analyzed for their VP 1 gene. 
All of the strains kept with the same cv or vv classification except one sample, Br/03/DR. This sample was 
classified as vv by its VP2 gene, but it was most closely related to the cv strains by its VP1 partial sequence 
and phylogeny. Studies on the phylogeny of VP1 have suggested a possible reassortment event that 
originated the vvVPl. In this case, the sample carrying vvVP2 and cvVPl could be a descendant of IBDV 
ancestors prior to the reassortment of wVPl; alternatively, it could be the result of a genetic exchange 
between the segments of different strains or with a live attenuated vaccine. Nevertheless, this is the first 
report of natural genetic reassortment of IBDV in Brazil. 

Key words: Brazilian infectious bursal disease virus; very virulent strain; VP1; sequencing; phylogenetic 
analysis. 



INTRODUCTION 

Infectious bursal disease virus (IBDV) causes a disease 
among young chickens and is of great economic importance to 
the poultry industry worldwide with regard to both mortality 
and immunosuppression (19, 20, 24). A member of the 
Birnaviridae family, IBDV contains a bisegmented dsRNA 
genome, designated A and B, within a non-enveloped 
icosahedral capsid of 60 nm in diameter. Segment A of the 



genome contains two partially overlapping open reading 
frames (ORFs). The first, small ORF encodes a 17 kD 
nonstructural host membrane -associated protein, VP5, and the 
second, large ORF encodes a 110 kD polyprotein that is auto- 
processed to produce the precursor of VP2 (pVP2), the 
multifunctional protein, VP3, and the viral protease, VP4. 
pVP2 is further processed to generate the outer capsid protein, 
VP2, which contains antigenic epitopes that are responsible for 
the induction of neutralizing antibodies. The segment B of the 
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genome encodes a 90 kD VP1 protein with RNA-dependent 
RNA polymerase activity (6, 16, 20, 24, 28). 

Two distinct serotypes of IBDV are recognized (1 and 2). 
Only serotype 1 is pathogenic for chickens and is classified 
according to the antigenicity and/or pathogenicity in classical 
virulent (cv) strains, very virulent (w) strains, antigenic variant 
strains, and attenuated strains (3, 25, 28). The antigenic 
variation has been based on mutations in a hydrophilic 
hypervariable region of the VP2 gene (4, 7, 20, 23). More 
specifically, only a few amino acids (aa) in this region 
contribute effectively to the antigenicity (28). Although many 
reports have also considered the aa of this same region of VP2 
as the determinants of virulence (1, 19, 20, 24, 28), further 
studies on the VP1 gene demonstrated the contribution of this 
gene to the virulence of IBDVs (2, 17). Both genes could be 
associated with the pathogenicity of IBDVs, yet the virulence 
markers remain unknown. Yamaguchi et al. (29) and Le Nouen 
et al. (14), studying the full genomic sequences from segments 
A and B, were able to define phylogenetically representative 
partial nucleotide sequences for use in molecular epidemiology 
studies. The purpose of this study was to analyze the partial 
sequences and phylogeny of Brazilian IBDV samples, as based 
on one of the regions of the VP1 gene, previously characterized 
as cv and w strains by the partial sequence of VP2 (8). 

MATERIALS AND METHODS 

Field Strains 

The bursae of Fabricius were collected from commercial 
layer and broiler chickens from Sao Paulo and Parana States, 
Brazil, between 1997 and 2004. The tissues were previously 
prepared, and identified as IBDV (8) (Table 1). 

Reverse transcription and polymerase chain reaction 
(RT/PCR) 

The viral RNA from the bursal homogenates was 
extracted using the QIAamp Viral RNA Kit (Qiagen Inc., 
Valencia, CA) according to the manufacturer's instructions. 



The primer pair used for the partial amplification of a 594 bp 
fragment of the VP1 gene (nts 290 to 883) is described in Le 
Nouen et al. (15). The amplification conditions were adapted 
from Hernandez et al. (10) for a one-tube procedure using the 
AccessQuick™ RT-PCR System (Promega, Madison, WI) and 
a Perkin-Elmer Co. (Norwalk, CT) thermal cycler. An aliquot 
of 5 \A of viral dsRNA was mixed with 4.5 jul of ultra-pure 
water, denatured by boiling for 5 min at 98°C, and immediately 
immersed in ice for 5 min. The following mixture was then 
added to the tube: 12.5 (j.1 of AccessQuick Master Mix (2x), 1 
jj.1 of each primer and 1 (j.1 of AMV reverse transcriptase. The 
cDNA was synthesized at 42°C for 60 min, and the 
amplification included an initial denaturation step (5 min at 
95°C) and 35 cycles of 95°C for 30 sec, 55°C for 1 min and 
72°C for 1 min. A final extension was performed at 72°C for 5 
min. The PCR products were subjected to electrophoresis 
through a 1.5% agarose gel supplemented with 1% ethidium 
bromide. 

Nucleotide sequence and phylogenetic analysis 

The PCR products were purified using a PCR Purification 
Kit (Promega) and were sequenced using Big Dye Terminator 
Cycle Sequencing Ready Reaction (Applied Biosystems) and 
an ABI 377 automatic DNA sequencer. The sequences were 
aligned using Clustal X and subjected to a BLAST search, and 
the phylogenetic analysis was performed using the distance 
based neighbor-joining (NJ) method with the Kimura two- 
parameter option in PAUP 4.0 with 1,000 bootstrap replicates. 
Nucleotide accession numbers of the samples and the strains 
used for the comparison are shown in Table 1 . 

RESULTS 

A fragment of the VP1 gene, between the clusters a and b, 
that is representative of the region used for phylogenetic 
analysis was successfully amplified in the seventeen Brazilian 
IBDV samples from our previous study (Table 1). The 
nucleotide (nt) and amino acid (aa) sequences were compared 
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among the samples and with those of the vaccine strains and other 
reported IBDV strains. The three previous cvVP2 IBDV samples 
were also classified as cv strains using the VP1 aa sequences. In 
thirteen of the samples identified as vvVP2, all of the four 
common and characteristic aa substitutions of the vv strains within 



the VP1 fragment were found (Table 2): the TDN tripeptide 
between aa 145-147 and 242E. The exception was the Br/03/DR 
sample that presented the conserved aa of the cv strains in these 
positions. 



Table L Description of infectious bursal disease virus strains used in this study . 



Strain 


Origin/Year 


VP1 Gene 








Pathotype 6 


Accession number 


Br/97/BA A 


Brazil/1997 


cv 


HQ602714 


Br/00/BS A 


Brazil/2000 


vv 


HQ602722 


Br/01 /BGK A 


Brazil/2001 


vv 


HQ602718 


Br/01 /BGS A 


Brazil/2001 


vv 


HQ602719 


Br/01 /BJM A 


Brazil/2001 


vv 


HQ602720 


Br/01 /BNK A 


Brazil/2001 


vv 


HQ602721 


Br/01 /CA A 


Brazil/2001 


vv 




Br/01 /CB A 


Brazil/2001 


vv 


HQ602723 


Br/01 /CC A 


Brazil/2001 


vv 


HQ602724 


Br/03/CK A 


Brazil/2003 


vv 


HQ602725 


Br/03/DB A 


Brazil/2003 




HQ602727 


Br/03/DQ A 


Brazil/2003 


vv 


HQ602728 


Br/03/DR A 


Brazil/2003 


reassortant 


HQ602717 


Br/03 /DU A 


Brazil/2003 


vv 


HQ602729 


Br/04/CR A 


Brazil/2004 


vv 


HQ602726 


Br/04/DE A 


Brazil/2004 


cv 


HQ602715 


Br/04/DN A 


Brazil/2004 


cv 


HQ602716 


UK-661 


UK/1994 


vv 


AJ878666 


OKYM 


Japan/ 1997 


vv 


D49707 


99009 


Brazil/99 


w 


AJ878672 


94432 


France/ 1994 


vv 


AJ878678 


variant A 


USA/1985 


variant 


AJ878675 


Cu-lwt 


German/1975 


cv 


AJ878644 


F52-70 


UK/1970 


cv 


AJ878643 


23/82 


Serotype 2/1985 


avirulent 


AJ878684 


02015.1 


Venezuela/2002 


reassortant 


AJ878685 


02015.2 


Venezuela/2002 


reassortant 


AJ878686 


ZJ2000 


China/2000 


reassortant 


DQ166818 


SH95 


China/2003 


reassortant 


AY134875 


02015.1 


Venezuela/2002 


reassortant 


AJ878685 


CEVAC IBDL 


vaccine 


attenuated 


AJ878660 


D78 


vaccine 


attenuated 


AJ878654 


228E 


vaccine 


attenuated 


AJ878657 


Bursine 2 


vaccine 


attenuated 


AJ878655 



^Brazilian samples obtained in this study. 
B cv: classic virulent; vv: very virulent. 



Table 2. Amino acid (aa) exchanges within the VP1 fragment (68-248aa) in different infectious bursal disease virus strains. 

aa at positions 



Strain 


Pathotype 


145 


146 


147 


242 


Cu-1 wt 


cv 


N 


E 


G 


D 


F-52/70 


cv 


N 


E 


G 


D 


Br/01 /BJM A 


cv 


N 


E 


G 


D 


99009 


vv 


T 


D 


N 


E 


UK661 


vv 


T 


D 


N 


E 


OKYM 


vv 


T 


D 


N 


E 


Br/03/DU A 


vv 


T 


D 


N 


E 


02015.1 


reassortant 


T 


E 


S 


D 


02015.2 


reassortant 


N 


E 


G 


T 


SH95 


reassortant 


N 


E 


G 


D 


ZJ2000 


reassortant 


T 


D 


N 


E 


Br/03/DR A 


reassortant 


N 


E 


G 


D 


variant A 


variant 


N 


E 


G 


D 


D78,228E and Bursine 2 


vaccine 


N 


E 


G 


D 


CEVAC 


vaccine 


N 


E 


D 


D 


23/82 


serotype 2 


N 


E 


G 


D 
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The identity percentages of the deduced aa sequence of 
the Br/03/DR sample with the others strains were determined. 
The sequence of this sample showed a higher identity with our 
Brazilian (Br) cv samples (93.6%-96.8%) than with the vv 
samples (91.0%-95.2%). However, the Br/03/DR sample had a 
91.0% similarity to the sequences of the cv reference strains, 
90.4%-91.0% to those of the vv strains and 91.5%-92.5% to the 
vaccine strains. Regarding the reassortant strains, our sample, 
Br/03/DR (vvVP2cvVPl), exhibited 92.0% similarity to the 



SH95 (wVP2cvVPl) and 90.4% to the ZJ2000 (cvVP2vvVPl) 
strains. The phylogenetic tree of the predicted aa sequences of 
the VP1 gene fragment from our samples and the published 
IBDV strains is shown in Figure 1. Our vv samples branched 
together with the other vv strains and in a distinct cluster from 
the cv strains, which included our cv samples and the 
Br/03/DR sample. This profile of the Br/03/DR sample 
carrying vvVP2, but cvVPl, is indicative of a genetic 
reassortment. 




Figure 1. Phylogenetic tree of 13 Brazilian IBDV isolates (rhombus symbol), reference strains (cv and vv) and vaccines. The 
neighbor-joining analysis was based in predicted amino acid sequences of the fragment of the VP1 gene (aa 68-248) using the 
PAUP program. The bootstrap values were included in the branches since that >50%. 
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DISCUSSION 

The fundamental molecular basis for the pathogenicity of 
IBDV is not entirely clear and remains questionable. To 
understand the circulating viruses in Brazil better, we 
continued our studies of the Brazilian samples from a period of 
seven years, which were previously classified as cv or vv 
strains by analysis of the VP2 gene (8), focusing on the VP1 
gene in this study. Through sequencing, comparative 
alignments and phylogeny of the VP1 gene, all samples 
maintained the previous cv and vv assignments, except one 
(Br/03/DR). 

Studies involving the multiple, full-length sequencing of 
the VP1 gene of diverse vv strains and from different countries 
have revealed the presence of 17 common aa residues, with 
eight that are unique to these strains (5, 10, 13, 18, 29). The 
Br/03/DR sample was the sole vvIBDV strain that did not 
present these common aa. Conversely, it presented the 
conserved aa of the cvIBDV strains, residue El 46, which Ren 
et al. (21) correlated to the attenuated virulence of the Chinese 
isolate, YS07, together with two other residues (41 and S687). 
However, these authors also correlated other unique mutations 
in the VP2 sequence to the low mortality of this isolate. The 
role of cvVPl in virulence of our sample cannot yet be fully 
established because of the absence of in vivo studies and full- 
length sequencing. 

The phylogenetic tree constructed from the partial 
deduced aa sequences of VP1 showed that the Brazilian vvVPl 
samples branched in a distinct cluster from the cvIBDV strains 
and even from the serotype 2 strain described in other studies 
(5, 12, 29). Phylogenetic studies with the VP1 gene of the 
vvIBDV strains suggested that the B segment is genetically 
distinct from that of all other IBDV strains, suggesting that the 
origin of the wIBDVs is a a genetic reassortment of an 
unidentified source (1, 12, 29). Using the estimation of the 
time of emergence of the most recent common ancestor 
(TMRCA) and coalescence techniques in both genome 



segments to the IBDV, Hon et al. (11) verified that this 
segment B reassortment occurred in the mid-1980s, whereas 
the emergence of the vv A segment was 20 years prior. The 
authors concluded that the enhancement of virulence and the 
expansion of the vvIBDV strains in mid-1980s may be due to a 
synergistic effect of mutations on both of its genome segments. 
According to this conclusion, the Br/03/DR sample presenting 
a B segment derived from a cvIBDV strain, but with the 
segment A of a vvIBDV strain (8), could be a descendant of 
the IBDV ancestors before the reassortment of segment B of 
vvVPl. 

It is known that the use of vaccines with intermediate 
virulence to control the disease by the vv strains can be 
responsible for a partial immunosuppression and mild bursal 
lesions. Although the viral multiplication does not revert the 
virulence, it could be the source of a genetic exchange 
producing a natural reassortant virus. Therefore, there is the 
risk of using live vaccines with the generation of such 
reassortments (26, 27). It could be another explanation for the 
origin of the Br/03/DR sample. Natural reassortments have 
mainly been described in China. The SH95 (20) and Gx (9) 
isolates and our Br/03/DR sample have vvVP2 and cvVPl. 
However, these two Chinese isolates maintained their virulence 
even without the vvVPl. Boot et al. (1) had initially verified 
that the VP1 gene would not have a major influence on the 
virulence of the vvIBDV strains. However, the ZJ2000 and 
TL2004 reassortants presented high virulence even without the 
presence of vvVP2 (cvVP2/wVPl) (26, 27). In this case, the 
VP1 gene could be contributing to the enhanced virulence of 
these wIBDVs, corroborating the further studies of Boot et al. 
(2). 

High mutation rates occur in RNA viruses, but these 
variations are limited by functional constraints, and the rate of 
genomic reassortment events in nature is also low and depends 
on many factors. Gao et al. (9), in a quasispecies population 
study of the Gx isolate, found a much higher frequency of 
mutation in VP1 than in VP2, an unexpected result for a RNA 
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polymerase. However, the most conserved region was 
correlated with another possible functional region. This isolate 
is also a reassortant, with VP2 derived from vvIBDV and VP1 
from a cv strain and a retention of virulence. Studying more 
than 50 strains of IBDV of diverse years and localities, Le 
Nouen et al. (15) concluded that the co-evolution of both the A 
and B segments of IBDV is frequent, although they have also 
found a natural reassortant (strain 02015.1). Continuing studies 
of full-genome sequences and the pathogenicity in vivo are 
very important to understand the true role of VP1 in virulence 
and/or of the risk of the use of live vaccines and that will 
provide the exact relationship between the genetic composition 
and pathogenicity. Regardless, this is the first report of a 
natural genetic reassortment of IBDV strains in Brazil. 
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